Analytic Solution of Ground State for 
Coulomb Plus Linear Potential 



W. Q. Chao(Zhao)i' 2, C. S. Ju^ 



1. China Center of Advanced Science and Technology (CCAST) 
(World Lab.), P.O. Box 8730, Beijing 100080, China 
2. Institute of High Energy Physics, Chinese Academy of Sciences, 
P. O. Box 918(4), Beijing 100039, China 



Abstract 

The newly developed single trajectory quadrature method is applied to solve 
the ground state quantum wave function for Coulomb plus linear potential. The 
general analytic expressions of the energy and wave function for the ground state 
are given. The convergence of the solution is also discussed. The method is applied 
to the ground state of the heavy quarkonium system. 



PACS: U.lO.Ef, 03.65. Ge 



Recently a new method has been developed by R. Friedberg, T. D. Lee and W. Q. 
Zhao [1,2] to solve the N- dimensional low-lying quantum wave functions of Schrodinger 
equation using quadratures along a single trajectory. Based on the expansion on 1/g, 
where g isa. scale factor expressing the strength of the potential, Schrodinger equation can 
be cast into a series of first order partial differential equations, which is further reduced 
to a series of integrable first order ordinary differential equations by single-trajectory 
quadratures. New perturbation series expansion is also derived based on this method, 
both for one-dimensional and N-dimensional cases. Some examples for one-dimensional 
problems have been illustrated in [1,2]. 

Coulomb plus linear potential has been widely applied to describe the heavy quarko- 
nium state. However, it is difficult to obtain an analytic expression of the energy and 
wave function. Here the single trajectory quadrature method is applied to solve the 
ground state for Coulomb plus linear potential. The general analytic expressions of the 
energy and wave function for the ground state are given. The convergence of the solution 
is also discussed. The result is applied to describe the ground state of heavy quarkonia. 
Some discussions of the limitation of its applicability is given at the end. 

Let us consider a unit mass particle moving in a central potential. The Schrodinger 
equation in the 3-dimensional space is expressed as 

[--V^ + V{r)]^!{r) = E^!{r). (1) 

Based on the single trajectory quadrature method[2], following steps should be taken to 

solve the problem for the ground state. 

1. For the potential V{r), the scale factor g is introduced as 

V{r)=g'v{r). (2) 

First, the energy is expanded in terms of 1/g: 

E = g'Eo + g'-'Ei+g'~^'E2 + --- (3) 

The highest g-power, I, in the expansion could be fixed by the dimensional consideration 
in the following way: Assume the behavior of the potential v{r) approaches when r —* 
and Eq is dimensionless. The dimension of in the first term of the right hand side in 
(|^) is the same as [r~^]. The dimension of each term in (|1|), namely [r~^] for V^, [g^r"'] 
for V{r) and [(?'] for E should be the same. This then gives 

. (4) 

n + 2 ^ ' 

The ground state wave function \E'(r) is expressed as 



2 



Substituting (|]) into (|l|) a equation for S{r) and E is obtained as 

^V^S{r)-^{VS{r)y + V{r)-E = 0. (6) 

Then S{r) is also expanded in terms of 1/g: 

S{r) = g^S,{r) + g^-^ S,{r) + g^'^^ S,{r) + ... (7) 

Substitute the expressions (H) and (j^) into the equation (||). By equating the coefficients 
of each g~^, a series of first order differential equations could be obtained. Now the highest 
(yf-power, m, in the S'-expansion should be determined. The highest power of g^"^ comes 
from the second term — |(VS')^ of the left hand side of (|^), which gives — |(VS'o(r))^ 
to the first one of the series of equations. The highest (^-powers in V{r) and E are g^ 
and g^ respectively. When k > I, v{r) should enter the first equation, as in the case of 
Harmonic Oscillator potential [1,2], which requires 2m = k; in the second equation for 
g^"i-j power, VS'o ■ VSi should be related to Eq, which gives 2m — j = I. On the other 
hand, when k < I, the first equation should be related to Eq and this gives 2m = /; the 
second equation with the power g'^^~i is then related to f (r), which gives 2m — j = k. 
From the above discussion we derive the following condition: 

2m = k, 2m ~ j = I, for I < k, 

2m = I, 2m-j = k, for I > k. (8) 

To ensure that {Ei} and {Si} enter the equations successively, we have i = j in (H) and 
(0). 

For the Coulomb plus linear potential 

V{r) = g\-- + fir) (9) 
r 

k = 2 and 

v(r) = h fir, (10) 

r 

which gives n = —1. From and (§) this gives I = 4 > k, which leads to m = 2 and 
i = j = 2. Therefore, for this potential (|^) and (|^ are expressed as 

E = g^Eo + g'E^ + E2 + g~'E, + --- + g-^''^-'^E^ + ---, 

S = g^So + S,+g-^S2 + g''S, + --- + g~^^^-^^S^ + ---. (11) 

Substituting (0) into equation (||), equating the coefficients of each g"' term, the following 
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series of equations are obtained: 

{VSoY = -2Eo, (12) 
VSo-VSi = ^V^So + v{r)-Ei, (13) 

VSo-VS2 = ^V25i-i(V5i)2-E2, (14) 
1 1 ""^ 

VSo-VSn = -V'S^_,--J2VSm-VSn-m-E^ (15) 

m=l 

2. Following ref.[2] the series of equations (|T^15) could be solved easily. Considering 
VSo = ^ the equation ([ID gives 



So{r) = ^-2Eo r. (16) 
Substituting (0) into (13) for Ei and Si(r), considering V^r = -, we have 



-2Eo^ = + (-- + /i r) - E,. (17) 

ar 2 ^ r r 



To keep Si{r) regular at r = we have 



This gives 



Eo = -l and So(r)=r. (1^ 



ar 

Si{r) = ^fir^-E,r. (19) 
Substituting (|^) into (14) for S'2(r) and E2, considering V^r^ = 6, we have 

'^ = lil^^■Q-'-^)-k^^r-E,r-E2. (20) 
ar 2 2 r 2 

In order that 5*2 (r) be regular at r = 0, we have 

E, = 0, 5i(r) = ^/i (21) 



This gives 



/iV^ + -/i - 



dr 2"^ 2 

1 o q 3 



2, 



^2(r) = -^^V + ^/i r - ^2 r. (22) 

D 2 



Following similar procedure for 5*3 (r) and E^, introducing V^r^ = 12 r, we have 

§ = l(-2,^ . + 5^1 - 2^) - ^ ri-l,V + I - B.) - E,. (23) 
In order to have S^^r) also regular at r = 0, we get 

^2(r) = -^A'- (24) 
o 



Then (Esf) becomes 



(i& o 1 



dr ^2 

53(r) = + - r. (25) 

Z o 



Continuing similar argument we could reach 



Es = 0, Ssir) = + 1 A', 

E. = S,ir) = -^/.V^ - (26) 



Now we introduce the general expression 



«{"V"-2V"-' (27) 



o<z<f 



dSfi 
dr 

then we have 

Sn = y 7 — — -«S"V"+^-2V"-' 

o<1^„ (n + 1 - 20 ^ 
V^5n = ;^r^§i = E ai-\n + 2 - 2/)r— V""' . (28) 

Substituting (^) and ( pSf ) into (15) and comparing the coefficients of the same power of 
r we obtain 

-^E E cf^atT-). (29) 
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To keep or S'2„,+i regular at r = we have 

E2n~l = 0, E2n = ^"S-Y'V"- (30) 

From o'q^ = 1 all the and En for ri > could be derived based on (pQ]) and (|30D. 
Combining with Eq = and = r we obtain 



1 1 1 

2 2 3 

2^^ +6?^ ^ +2?' 



W(r) = exp[-5( r - -fi r + — r + — r 

1 Q /I 3 o n 1 



The same procedure can be performed if we define 

e = g'^n (33) 
and solve the equation (|ip for the potential 

V(r) = -— + er. (34) 
r 

The derivation has been given in ref.[2] and the result is exactly the same as (^) and 
(32), considering (^) we have: 



iTr/N r2 -'- 2,-^23,-^22 

^(r) = exp[-g r ~ —t r + — e r + — e r 

1 o ^ 3 o o , 1 



,,eV -e^ r'^ + — -e^ + ■ ■ ■] . (36) 

8^10 4^12 8^14 J ^ I 

Introducing a parameter X = e/g^ = fi/g^ the energy could be expressed as 

i^ = ^?Vi + E?"iVA"). (37) 



2 ^2 

n>l 



Introducing e„ = ^ctn-i the energy could be expressed as 

E = g\~ + Y.enXn. (38) 

n>l 

This method can easily calculate the energy expansion series up to any order of n. It 
gives the possibility to analyze the convergence of the series in details. The convergence 



6 



of the expansion series of the energy E depends on the parameter A. In fact, this series is 
an asymptotic one. For certain value of A we could only reach a certain accuracy of the 
energy. In Table I the ratio of i?„ = \en/en-i\ for different n is hsted. From (|37|) we know 



that the ratio of the successive terms in the energy expansion series is A|en/e„_i| = XRn- 
It can be seen that the series would be meaningful only when A < This gives the 
limitation of the applicability of this method. It also tells us how accurate the final result 
could reach for a fixed value of A. In table II for each order n, the A corresponding to 
XRn ~ 1 is listed. The obtained energy E/g^ for this special value of A at each order 
of n is also given, together with the reached accuracy e„A". For example, when n = 11 
the corresponding A = .052 gives Ai?ii ~ 1. It means that the correction term increases 
when n > 11 and increases further. This would finally give a divergent result. Up to 
n = 11 the obtained CnA" = 0.00004 which gives the accuracy of the obtained energy 
E/g^ = -0.43443 for A = 0.052. 

Now consider a pair of heavy quark and antiquark with equal mass m and color charge q 
and —q, moving in a Coulomb plus linear potential. In the non-relativistic approximation 
the wave function ?/'(r) to describe the relative motion of the quark pair satisfies the 
following Schrodinger equation: 

1 0^ 

r -V^- — + Kr]V^(r) =^V^(r), (39) 



2 • m/2 r 

where k is the strength of the linear potential. After a simple transformation, (|39D becomes 
(0) and with 

-^mq = g , 

-mS = E, 
2 

1 1 
— -rriK = u or —rriK = e 
2g^ ^ 2 

^ = —n^ = ^ = 4- 40 

'm Q 9 9 



Substituting (^) into (|37D the ground state energy £ of the two quark system can be 
expressed as 

- \m<,H-\4x-\x'.?lx^...). (41) 

For the quarkonium system the color charge q^ can be related to the strong coupling 
constant as [3] 

= \as. (42) 
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For example, taking = 0.5 and n = 1 GeV/fm we have 



A = ^ ~ 6AGeV^ x At- (43) 



Taking the mass of the charm quark rric =1.6 GeV, we have 



Ac = ~ 2.5 , (44) 



which is too large for the application of this method. If we look at the expansion series, 
for the ground state of J/ ip we have 



mj/^ = 1.6 Gey X 2 + 

~ 3.2 Ge1/ + .2 Ge\/(-.5 + 3.7- 9.4 + 105 + ---). (45) 

Obviously, this expansion has a very bad behaviour. Therefore this method could not be 
applied to charmnium since the charm quark mass is not heavy enough. For bottom and 
top quark, taking mb ~ 4.5 GeV and rrit ~ 170 GeV , we could obtain 

Aj, = 4^ ~ 0.32 and At = ~ 2.2 x 10"^ (46) 

For the ground state of the bottom and top quarkonium we have 

m^l = 4.5 GeV x 2 + 

~ 9.0 Ge\/ + . 56 GeF(-.5 + 0.48 -0.15 + 0.22 + ■■■) (47) 

rritj = 170 GeV x 2 + £b 

~ 340 Ge1/ + 21.3GeV(-.5 + 3.3 X 10"^ -7.3 X 10"^ + ■■■). (48) 

From d^ ) and (48) it can be seen that in the case of bottom quarkonium A;, ~ 0.32, the 
accuracy of the derived energy has been improved and this method could well be applied 
to the top quarkonium. 
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Table I. Rn = |en/en-i 



n 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


Rn 


4.5 


7.36 


9.67 


11.62 


13.35 


14.95 


16.45 


17.90 


19.32 


20.72 




.22 


.11 


.10 


.086 


.075 


.067 


.061 


.056 


.052 


.018 



Table II. The obtained energy and its accuracy for different A 



n 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


A 


.22 


.14 


.10 


.086 


.075 


.067 


.061 


.056 


.052 


.048 


E/g' 


.20 


.31 


.361 


.379 


.394 


.4048 


.4130 


.4198 


.42535 


.43089 


|e„|A" 


.07 


.02 


.006 


.002 


.001 


.0004 


.0002 


.0001 


.00004 


.00002 
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